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Abstract The asymmetric exclusion process (ASEP) has attracted a lot of interest not 
only because its many applications, e.g. in the context of the kinetics of biopolymerization 
and traffic flow theory, but also because it is a paradigmatic model for nonequilibrium sys- 
tems. Here we study the ASEP for different types of updates, namely random-sequential, 
sequential, sublattice-parallel and parallel. In order to compare the effects of the differ- 
ent update procedures on the properties of the stationary state, we use large-scale Monte 
Carlo simulations and analytical methods, especially the so-called matrix-product Ansatz 
(MPA). We present in detail the exact solution for the model with sublattice-parallel and 
sequential updates using the MPA. For the case of parallel update, which is important 
for applications like traffic flow theory, we determine the phase diagram, the current, and 
density profiles based on Monte Carlo simulations. We furthermore suggest a MPA for 
that case and derive the corresponding matrix algebra. 
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1 Introduction 



For nonequilibrium systems in low dimensions an understanding can often be gained by 
studying rather simple models []!]-§). The one- dimensional asymmetric exclusion process 
(ASEP) has been used to describe various problems in different fields of interest, such as 
the kinetics of biopolymerization |7| and traffic ||. On the other hand, the ASEP is so 
simple that it has achieved a paradigmatic status for nonequilibrium systems ||. It can 
be mapped onto a surface growth model known as the single-step model [TtJ and in the 
appropriate hydrodynamical limit its density profile obeys the Burgers equation which is 
itself closely related to the KPZ equation [fLlp . The ASEP can also be viewed as a prototype 
for so-called boundary- induced phase transitions fL2f: the boundaries, which can inject and 



remove particles from the system, govern - in a subtle interplay with the local dynamical 
rules - the macroscopic behaviour of the model and can produce different phases and phase 
transitions. 

In |13| recursion relations on the system size have been derived for the ASEP with random- 
sequential update and open boundary conditions. Open boundaries here and in the fol- 
lowing mean that particles are injected at one end of a chain of L sites with probability a 
and are removed at the other end with probability (3. In [|T3|] these recursion relations were 
solved for the special case a = (3 = 1. This solution was later extended in |14j] to general a 
and f3. A very elegant solution of the general case was given at the same time in [15] using 



a matrix product Ansatz (MPA) for the weights of the stationary configurations similar 
to the matrix product groundstate of certain quantum spin chains ||16|| . This Ansatz can 
be used to compute density profiles and correlation functions. The relationship between 
the MPA for quantum spin chains and one-dimensional models of statistical physics will 
be discussed in this paper. 

Since then the MPA was extended to find the transient of the model 



171, to describe 



the ASEP with a defect in form of an additional particle with a different hopping rate 
|T8|1 or a blockage fl9fl , to solve the case of oppositely charged particles (with hard-core 
repulsion), which move in opposite directions (driven by an external electric field) and can 
interchange their charge if they meet p0|. The MPA was also used to recover solutions of 



certain integrable reaction-diffusion models |2T 



Most of these solutions have been found for random-sequential dynamics. In that case the 
master equation can be rewritten as a Schrodinger-like equation for a "Hamiltonian" with 
interactions between nearest-neighbours only [22]. Other updates lead to more complicated 
master equations with non-local interactions. 

It also turned out that the MPA is not just an Ansatz. The stationary state of an one- 
dimensional stochastic model with arbitrary nearest-neighbor interactions and random- 
sequential update can always be written as a matrix product p^| . This is also true for 
an ordered-sequential or sublattice-parallel update, which can be shown to be intimately 
related ||24|| . The matrices in the MPA are generally infinite-dimensional. Therefore, eval- 
uating physical quantities such as density profiles is still a formidable task, but often it 
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is at least possible to obtain asymptotic expressions. However, in certain regions of the 
parameter space the matrices can reduce to low- dimensional variables. A simple example 
is the case of one-dimensional matrices, which is equivalent to a mean-field solution. 

The implementation of the order of application of the local transition rates for a given model 
(the type of update) is an essential part of the definition of the model, since the transient 
and even the stationary state may differ dramatically ||25|| . In the following, the ASEP will 



be studied with types of update that are often more useful for Monte Carlo simulations 
than the random-sequential update. The aim of this paper is to use the ASEP as a case 
study to investigate the consequences of different types of updates onto the stationary state 
of a nonequilibrium model with open boundaries. Concerning the analytical treatment, the 
MPA turns out to be a useful tool. 



Hinrichsen [26] was the first to apply the MPA to the ASEP with sublattice-parallel update 
and deterministic bulk dynamics. He could confirm earlier conjectures for the correlation 
functions f27j. It is important to note that this update is substantially different from 



the fully-parallel update used for example for modeling traffic flow. Nevertheless, to our 
knowledge this is the first model with a discrete-time update which has been solved using 



the MPA. The results presented here build partially on our previous work ||28[| , where we 
found a mapping of the ASEP with ordered-sequential and sublattice-parallel update onto 
the random- sequential case. For sublattice-parallel update this was done independently by 



Honecker and Peschel [J2I| . We will explain the mapping in more detail and give a physical 
interpretation of the underlying Ansatz. As a by-product, we solve the model on a ring. 

We also study numerically and analytically the ASEP with parallel update and present new 
results like the phase diagram and the current. We show that the model can be described by 
a matrix-product structure at least for small systems. The underlying algebra is somewhat 
different to the other ones because a third state appears, which is neccessary to decompose 
the transfer matrix into a more simple product. 

This paper is organized as follows: Section |^ gives the definitions of the model and the 
updates used and establishes a precise relationship between the MPA for the stationary 
state of a stochastic model and for the groundstate of quantum spin chains. Section || 



is devoted to the MPA for the ASEP with discrete-time updates. In Section |3J] the 
matrix algebra for the ordered-sequential and the sublattice-parallel update is derived. 
In Section |3.2| we propose a MPA for the parallel update by mapping the model onto a 



3-state model with an ordered-sequential update. Furthermore, we find a special line in 
parameter space where the 2-cluster approximation becomes exact. In Section [|we use the 
mathematical results of Section |3] and Monte Carlo simulations to calculate the current, 
density profiles and correlation functions for the discrete-time updates. Finally, Section |5| 
contains a concluding discussion. We have included several appendices, which mainly list 
details of the calculations. 



2 Definition of the ASEP and the MPA 
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2.1 The ASEP and the different updates 



Consider sites located on a chain of length L. Each site % (1 < i < L) is either occupied by 
a particle (r, = 1), or it is empty (rj = 0). We have three basic mechanisms (see Fig. 1): 
Hopping: We look at the pair (i,i + 1). If we find a particle at site i and no particle 
("hole") at site i + 1, we move the particle one site to the right with probability p. In the 
remaining cases, nothing happens. 

Injection and removal: At the boundaries particles/holes can be injected and/or removed. 
At the left boundary a particle can be injected with probability a if site 1 is empty. At 
the right end of the chain a particle will be removed with probability (3. 
The model can be generalized to allow hopping in both directions by introducing a prob- 
ability q for hopping onto an unoccupied site to the left. Furthermore, one can also inject 
particles at the right end and remove particles at the left end. These modifications do not 
change the basic features of the model and will not be considered in the following. 

We now have to define the order in which to perform the hopping, the injection and the 
removal in terms of time and space. There are four basic ways to do that: 

(a) random- sequential update: We pick at random a site i. If 2 < i < L — 1, each particle 
has a probability pdt of jumping to the right (if this site is empty). If i — 1 , we allow for 
hopping to the right with probability pdt and particle injection with rate a if this site is 
empty. For i = L a particle is removed with rate (3 if the site is occupied. This update 
is the realization of the usual master equation in continuous time. A different p would 
simply result in a rescaling of time (see Sect. As a consequence, the phase diagram of 
the ASEP with update (a) depends only trivially on p. Therefore one can set p = 1 which 
is most efficient for computer simulations. 

The following three updates are discrete in time. 

(b) sublattice-parallel update: We first use our rule for injection (removal) at site 1(L). We 
then perform our rules for hopping on the pairs (2, 3), (4, 5) etc. After that, we update the 
pairs (1,2), (3, 4) etc (L has to be even). This update can be used efficiently for computer 
simulations. Its main advantage for theoretical purposes is that its transfer matrix can be 
written as a product of local terms. 

(c) ordered- sequential update: We start at the right end of the chain and remove a particle 
at site i = L with probability j3. We then update the pair (i — L — 1, i — L). We continue 
with pair (i = L — 2, i = L — 1) and so forth, until the left end of the chain is reached. After 
the update of pair (i — l,i — 2), we allow for injection at site i — 1. For models where 
particles only hop to the right this update may also be called more precisely 'backward- 
ordered- sequential update'. 

Obviously, the order of update can also be reversed. For the ASEP, these two updates are 
connected by a particle-hole symmetry: injecting particles can be regarded as removing 
holes, and vice versa. Therefore it is sufficient to study just one of the ordered-sequential 
updates. 

We would like to point out that in principle one has to distinguish two different types 
of ordered-sequential updates which one could name site-ordered-sequential and particle- 
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ordered-sequential, respectively. In contrast to the site-ordered-sequential update described 
in (c) above, in the particle-ordered-sequential update the rules are only applied to occupied 
sites, i.e. to particles. This might have a strong effect, as can be seen most easily for the 
case of small particle numbers and the completely asymmetric deterministic case p — 1. In 
the case of site-sequential update a single particle injected at the left end moves through 
the lattice in one timestep (i.e., one sweep through the lattice). For the particle-sequential 
update a timestep means updating occupied sites only and so the particle moves only one 
site. By looking at a lattice with two particles, one can already see that the two different 
updates might introduce rather different correlations. Starting with particles separated by 
d empty sites, in the site-ordered-sequential update the left particle will move to the right 
until it reaches the right particle, which then starts to move. On the other hand, in the 
case of particle-ordered-sequential update the particles will stay always d or d — 1 sites 
apart. For general values of p the situation is similar. 

The differences between these two sequential update procedures manifest itself also in the 



solution for periodic boundary conditions. We will come back to this point in Section 4.1 



In the following we will always consider the site-ordered-sequential update — until stated 
otherwise — to which we will refer as ordered-sequential update for brevity, 
(d) parallel update: The rules for hopping, injection, and removal are applied simultaneously 
to all sites of the whole chainF]. 

The parallel update usually produces the strongest correlations and is used for traffic 
simulations ||. In the case of the ASEP, it is nearly identical to the particle-ordered- 
sequential update. 

Fig. 2 illustrates the updates (b) and (c). 

If the ASEP with updates (a)-(c) is put on a ring (no injection/removing of particles and 
periodic boundary conditions), a trivial stationary state (where correlations are absent, see 
below) is reached, while update (d) produces a particle-hole attraction H36fl . This already 
shows that different updates might yield a different behaviour. 

For analytic calculations usually the random-sequential update is most convenient since 
it can be formulated in terms of a "Hamiltonian" with nearest-neighbour interaction. In 
Monte Carlo simulations, however, ordered-sequential updates can be implemented more 
effectively. 

At this point it is necessary to point out the existence of some confusion in the nomenclature 
of the different updates. The random-sequential update (a) is sometimes simply called 
'sequential update', as is the ordered-sequential one. In several publications the sublattice- 
parallel update (b) is called 'parallel' which makes it necessary to refer to the parallel 
update (d) as 'fully-parallel'. We urge the reader to check carefully which type of update 
is actually used when consulting the literature. In order to avoid further confusion we will 
use the terminology which seems to us the most precise. 



1 Note that for the parallel update the partially symmetric ASEP (hopping to right and left with 
probabilities p and q, respectively) might lead to ambiguities. One therefore has to set q = 0. 
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2.2 Master equation and quantum formalism 

Our starting point is the master equation for an arbitrary one- dimensional stochastic pro- 
cess. Following we rewrite this equation as a Schrodinger-like equation in imaginary 
time. We consider a chain of L sites j with state variables Tj = 0, 1. Generalizations to 
the case where the state variables can take more than two values are straightforward. A 
configuration of the whole system will be denoted by {r} = {ti,t 2 , ...,tl}, its weight by 
P({rM). 

The master equation then has the form 

J^({t}, f) = £ [w(r' - r)P({r'}, t) - w(r -> r')P({r}, f)] (1) 

where w(r — > r') denotes the rate for a transition from {r} to {r'}. 

Let us now define a ket state \P(t)) in the following way. We take an orthonormal basis 
in the configuration space {t}, 

It) = fa, ....ti) (2) 



with (t'\t) = 5 T ' tT and define 



P(t)) = 5>({T},t)|T), 

{r} 



{t\H\t') = -w(t' -> r) 



(r|W|r)= £ w(r->r') 

{r}^{r'} 



(3) 

(4) 
(5) 
(6) 



i.e. we have P({r},t) = (r\P(t)). It is then easy to see that 

^\p(t)) = -n\p(t)) 

holds, with the (generally non-hermitean) "Hamiltonian" 7i given by 



for the off-diagonal elements (r ^ t') and by 



for the diagonal elements. 

From (|j) one can see that the stationary state \Pq) of the stochastic model corresponds to 
the "groundstate" with "groundstate energy" zero (Eq = 0) of the "quantum spin chain" 
defined by (|), (§), i.e. 

n\P ) = 0. (7) 
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The conditions (0), (^) guarantee that the real parts of the other eigenvalues E\ of TL 
are non-negative. Using the bra groundstate (0| = J2{ T }( T \i the avera g e of an observable 
*4.({t}) at time t is given by 

(A)(t) = £>({r})P({r},t) = (0\A\P(t)}. (8) 

M 

Expanding the initial state \P(t = 0)) = J2\a\\?P\) in terms of the eigenkets I^a) with 
eigenvalues E\ of 7i this can be rewritten as 



(A)(t) = (0\A\P(t)) = (0\Ae- nt \P(t = 0)) = J2 



a x e 



-E x t 



(o\AM. 



(9) 



This shows that the behaviour for large times t is governed by the low-lying excitations. 
We now restrict ourselves to stochastic processes with random-sequential dynamics and 

t't' 

local and homogenous transition rates r r J T ^J (denoting the rate for a local transition of 



sites j and j + 1 from {t'-t'- + ^) to 



)), which do not depend on time. In this case Ji 



is a sum of local Hamiltonians with nearest-neighbour interaction only |22 



(10) 



with 
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plO 
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1 00 
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1 01 

_ru 

1 10 

r ii _j_ r oi , pii ; 

1 00 + 1 01 ^ 1 10/ 



(the basis is (00), (01), (10), (11)). Each column of hj adds up to 0, because the probability 
has to be conserved. 

For the updates (b)-(d) it is generally not possible to write Ti, in the form flTU| ) since one 
would have to include terms acting on sites which are not neighboured. It is then more 
suitable to use directly the transfer matrix T describing the update of the whole chain 
during one timestep. In the case of the ordered- sequential update, T is simply a product 
of the local transfer matrices = hi — 1. It follows 



\P(t + l))=T\P(t)), 



(12) 



which means that in order to find the stationary state, one has to solve for the eigenvector 
of T with eigenvalue 1. 



2.3 Optimum groundstates and MPA for quantum spin chains 

The construction of optimum ground states for quantum spin chains via matrix products 
was introduced in [T(| (see also |§0|,^T[ and for further references). Let us consider 
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a Hamiltonian for a quantum spin chain (with periodic boundary conditions) of the form 
H! = Y^j=i hp where h'j is the local hermitian Hamiltonian and independent of j, acting 
only on spin j and j + 1. 

It is always possible to set the lowest eigenvalue of h'j equal to zero by adding a suitable 
constant. Then h'j is positive-semidefmite and since Tri! is the sum of positive-semidefmite 
operators, it follows that zero is a lower bound for the ground state energy Eq of Ti,', i.e. 
Eq > 0. Usually, E is greater than zero (E > 0) and the global groundstate involves also 
excited states of h'j. Therefore, a construction of the global groundstate is usually very 
difficult. 

However, there are special cases where E is equal to zero, 

ft'lVo) = (13) 

and therefore for all j 

K |-0o> = 0. (14) 



A state |-0q) is called optimum groundstate of TC if and only if condition ([14]) holds. This 
implies that the groundstate energy is independent of the system size, i.e. there are no 
finite-size corrections. 

The idea is now to construct groundstates by means of a product of matrices, 

\i>o) = Tr (mi <g> m 2 <g> ... <E> m L ) , (15) 

where the entries of matrix rrij are spin-1 single-site states and the symbol <g> denotes 
the usual matrix multiplication of matrices with a tensor product of the matrix elements. 
Note that rri\ (g> m 2 <8> ... <8> ttil is still of the same size as the matrices rrij, but its elements 
are large linear combinations of tensor product states. The trace assures the translation 
invariance of the groundstate. For non-periodic boundary conditions it has to be replaced 
by a suitable linear combination of the elements of m\ ® ® ... ® m^. 
As an example |16 |, the ground state of a large class of antiferomagnetic spin-1 chains can 



be constructed using the S z eigenstates |0)j and |±)j by 



m, - ( a ®> 'tl' ) , (16) 



-h d\a) 



where the a, b, c, d are real numbers. Condition ([14]) requires 

h'j {rrij ® m j+ i) = , (17) 
i.e. all four elements of rrij <g> rrij+i are local groundstates of hj. Let us now write 

m j =A o -\0)+A--\-)+A+-\+) (18) 
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with suitable 2x2 matrices A and A±. The '•' denotes a product of each entry of the 
matrix to the left with the single-site state on the right. 
It is obvious that (O) can be written equivalently as 



|Vo) = Tr 



or, defining a column vector 



A v 


/ A 




5 A. 


A + ' 





... <g> 



A 



(19) 



.4 



as 



|V>o) = Tr(v4 ® i ® ... A) 
The condition (|i~7|) then be rewritten as 

^. (A <g> A) = . 



(20) 



(21) 



(22) 



This means that there are two equivalent ways to write |^o)- While ( |T5D uses a product of 
matrices with vectors as entries (the usual notation for quantum spin chains), (|2lf) expresses 
|^o) as a product of vectors with matrices as entries. The original idea of Derrida et.al. 
15] was to construct the stationary state \P ) of a stochastic process defined by (|l(J) as a 



suitable linear form of a product of matrices where each matrix corresponds to a single-site 
state precisely as in (pTF). 



2.4 MPA for stochastic systems 

In the following we want to describe how the MPA can be applied to stationary states \Pq) 
of stochastic systems. In contrast to the case of quantum spin chains we already know 
the corresponding "groundstate energy" of the stochastic Hamiltonian Ti defined by (§) 
and (|5|). From (|) we see that 7Y|P ) — and therefore the groundstate energy is zero, 
independent of the system size. This hints at the applicability of the MPA and a possible 
generalization of the optimum groundstate concept. 

Indeed it turns out that, using an MPA of the form (0), (|i~4l) can be replaced by a more 



general condition by allowing for a divergence-like term |21[|, i.e. 



hj(A ® A) = A® A - A® A , (23) 

where the vector A can be different from the vector A. It is easy to see that (for periodic 
boundary conditions) the divergence- like terms cancel after summing over j. Hence, the 
stationary state of the stochastic process described by (^) and (^|) is of the form (f?T|). 
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For hermitian Hamiltonians one always has A = A and reduces to Therefore one 
can regard (^) as the generalization of the optimum groundstate concept to non-hermitian 
Hamiltonians. 

Again it is possible to generalize these ideas to treat non-periodic boundary conditions. As 
an example we briefly review the solution of the ASEP with random-sequential update and 
open boundary conditions. The "Hamiltonian" reads (hopping rate p, feeding and removal 
rates 

L-l 

H = h 1 + h L + ^h j (24) 



with the boundary terms 



a 0/ ' n \0 P 



and the bulk "Hamiltonian" 



(0 








o\ 








-p 











p 





Vo 








0/ 



(26) 



From (0) we can see that p only rescales time (and a and (3) and that it would therefore 
be sufficient to study p = 1 . 

Following the MPA, the matrix for the particle (hole) is denoted by D (E), so that A = 
E \ 

Since the ASEP with open boundary conditions has (in general) no translational 



D 

invariance, the trace in ( |21]) is replaced with a scalar product: 

\Po) = -^-((W\A®---®A\V)), (27) 

where the normalization constant is equal to Zi, = (W\C L \V) with C = E + D. 
The brackets ((...)) indicate that the scalar product is taken in each entry of the vector 
A®...® A. 

More explicitly (|27|) means that the weight P(n, . . . , Tl) of a configuration {r} in the 
stationary state is given by 

1 L 

P(n, . . . ,r L ) = —{W\ J] [r 3 D + (1 - rj )E] \V). (28) 

Thus one has a simple recipe for the calculation of the (unnormalized) weight of an arbitrary 
configuration {r}: Translate the configuration into a product of matrices by identifying 



10 



an empty site (tj = 0) with E and an occupied site (tj = 1) with D. For example, the 
configuration 011001 • • • corresponds to the product EDDEED ■ ■ ■ = ED 2 E 2 D ■ ■ ■ . Then 
multiply by the vectors (w\, \v) from the left and right, respectively 

We assume that A in (^3| ) is of the form A = y q j where E and D denote matrices 

acting in the same vectorspace as E and D. Executing the sum ([[(]) in (p3| ) leads to a 
cancelation of all terms in the bulk of the chain. The remaining terms at the boundaries 
vanish if the vectors (W\ and \V) are chosen appropriately: 

(W\h 1 A=-(W\A , h L A\V)=A\V). (29) 

Inserting (ESI) and (EBf) into (1231), we get a system of quadratic equations in E,D,E 



and D which is called the algebra of Ansatz (p3|). The dimension of the matrices is not 
determined by the Ansatz. However, taking E = 1 = —D, the equations reduce to the 
"DEHP algebra" |T§ 

pDE = D + E , (30) 
a(W\E = (W\ , (31) 
PD\V) = \V) . (32) 



Equation fl3"0|) can be guessed intuitively: the current Jl(j) (describing the flux of particles 



through bond j of a chain of length L) has to be constant throughout the chain. Jl(J) is 



given by JT5] 



JM = ■^-(W\C'- 1 DEC L -'- 1 \V) , (33) 

and we see that DE oc C = E + D is the simplest way to achieve a constant current. 

It is possible to derive explicit expressions for D, E, (W\ and \V) It turns out, that 

one finds one-dimensional representations if and only if 

a + (3 = p. (34) 

In all the other cases, the matrices are infinite-dimensional. 

Up to this point, the MPA appears to be just an Ansatz for the stationary state. However, 
it can be shown (2^| that the MPA ( 27|) with the "cancellation-mechanism" ( p3|) is an 



equivalent reformulation of the master equation for a stochastic process with random- 
sequential update and nearest-neighbour interaction. Therefore it is of general interest 
to study the quadratic algebras which are produced by (|23| ) and to try to find explicit 
representations ||33|-|35|. 

In the next section, we will construct the stationary state for the ASEP with sequential 



and sublattice parallel update. Note that we do not have a "Hamiltonian" of the form (10) 
in this case and that the "cancellation mechanism" will not be appropriate. 
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3 MPA for the ASEP with updates in discrete time 



In this section we will generalize the results for the random-sequential update to discrete- 
time updates. In Section |3.1| we solve the ASEP for sublattice-parallel and sequential 
updates using the MPA. In Section ^]2| we conjecture a MPA for the parallel update. 

Let us briefly discuss the precise connection between the random-sequential update and the 
ordered-sequential update, say from the right to the left, for an arbitray one-dimensional 
stochastic process with nearest neighbor interaction. We have seen, that in this case 
the " Hamiltonian" describing the random-sequential dynamics is of the form (10), TC = 
^2j=i hj, while for the ordered-sequential udpate one has to use the transfer matrix 



n 



(35) 



with local matrices tj connecting the sites j and j + 1. From (|T0D and (|9]) it is clear that 
one of the parameters can be used to rescale the time unit; in the case of the ASEP this 
leads to expressions which are independent of the hopping probability p. For the update 
(j35|) with discrete time, this is obviously impossible and we expect the phase diagram of 
the ASEP to depend on p nontrivially. Only in the limit of vanishing densities these two 
updates can be mapped onto each other |2{J by inserting tj — hj — 1 (where 1 denotes 



the identity matrix) into ( |3"o] ) and expanding the product. For nonvanishing densities, the 
equations (|T0D and (|3~5p are connected nontrivially via additional, non-locals terms. 



3.1 Sublattice-parallel and ordered-sequential update 

We now solve the ASEP with ordered-sequential (update (c)) and sublattice-parallel (up- 
date (b)) dynamics. So far, the ASEP with the latter update and deterministic hopping 
has been studied by Schiitz |27] and Hinrichsen [ 2f| . A brief account of our work has been 



given in p8 |. For simplicity, we will first concentrate on the ordered-sequential update (c) 
from the right to the left. Since this update is discrete in time, a stochastic "Hamiltonian" 
of the form (|10D is not at hand. Therefore, the transfer matrix T<_ has to be used. This 



means by definition 

|P(t+l)> = T^\P(t)) (36) 

The stationary state \Pq) must not change under the action of the and therefore is 
eigenvector with eigenvalue 1 of T<_: 

|Po)=T_|P ). (37) 

Let us now explicitly write down The boundary conditions can be represented by 
operators 1Z and C acting on site j = L and j = 1, respectively: 

M 1 ;"!)- (38) 
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The basis chosen for 1Z and £ is (0, 1). The update-rule for any pair of sites + 1) can 
be written as 
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(39) 



(40) 



L = £®1®...® 1, 
R = 1 <g> ... ® 1 ®TZ, 
Tj = 1 ® 1... ® 7} ® 1. 



(41) 



where 1 denotes the identity matrix. 

Formally, a "Hamiltonian" can be defined as = 1 — T«_. However, cannot be 
written simply as a sum of local "Hamiltonians" . This means that while we can try a MPG 
Ansatz fl2?D we cannot use the "cancellation" -mechanism (|23|). 
However, the sequential nature of T<_ suggests another mechanism: 



T 



A® A 



A® A. 



(42) 



(W\CA = (W\A, TZA\V) = A\V), (43) 

where [ J w ^ square matrices E,D. 

This means that a "defect" A - corresponding to a local perturbation of the stationary 
state defined by (p5TJ) - is created in the beginning of an update at site j = L, which is 
then transported through the chain, until it reaches the left end and disappears. 
Equation (|42[) leads to the following bulk algebra: 

[E,E] = [D,D] = 0, 

ED + pDE = ED, (44) 

(l-p)DE = DE, 

and the boundary conditions 

(W\E(l-a) = (W\E, (1-P)D\V) = D\V), 
(W\(aE + D) = (W\D, (E + (3D)\V) = E\V). 
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The ordered-sequential update in the opposite direction (left to right) can be treated in 
the same way. The stationary state is given by 

\P )-> = ^-((W\A®---®A\V)) (46) 



with the same mechanism (fE|) and fl43|) . However, it is more convenient to use the particle- 
hole symmetry for the calculation of averages. 

For completeness let us briefly discuss the sublattice-parallel update (b). In this case the 



transfer matrix has the structure 26. 2£ 



T s _|| = T x ■ T 2 (47) 

with 

Ti = c®t ®t • • • ® r ®n, 

T 2 = T ®T ■■■®T ®T, (48) 
and the MPA for the stationary state is of the form 

|Po)hi = -tt((W\A® A® A® A---® A® A\V)). (49) 

One can now use exactly the same mechanism (f42|), (f43|) as in the ordered-sequential case 



and thus obtains the same algebra (f44l) , (45). For p = 1 the algebra was first derived and 



solved by Hinrichsen [2B| using a two-dimensional representation for E, D, E, D, (W\, 
\V)- 

The fact that the ordered-sequential and the sublattice-parallel update lead to the same 
algebra flUD , (|45|) implies the existence of an intimate relationship between the averages of 
observables. Although the stationary states themselves are different, they are connected 
via transformations, and it can be shown that the density profile of the ordered-sequential 
update from the left (right) to the right (left) corresponds to the density of the even 
(odd) sites produced by the sublattice-parallel update |24j . This result holds for arbitrary 
stochastic models with nearest neighbor interactions. 
One can check that for 

(l- a )(l-/3) = l_p , (50) 

a one- dimensional solution of the algebra exists^]. This equation defines a line in the phase 
diagram where the mean- field solution becomes exact. It turns out that this line touches 
all phases. This makes it possible to calculate quantitites such as the current very easily, 
because the analytic expression for the current does not change inside a phase. 

For general values of a, (3 and p the algebra fl4"4] ) and ( fl"5D can be mapped onto the 
generalized DEHP-algebra [28|. This is shown in App. |A| and will be used later in Section 



App. [FJ discusses the consequences of the particle-hole symmetry for the two ordered- 
sequential updates in more detail. App. deals with the special case of symmetric diffusion. 

2 If wc allow hopping in both directions, this line is given by |28j (1 — a)(l — /3)(1 — q) = 1 — p where q 
is the hopping probability to an empty site on the left. 
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3.2 Parallel update 



As far as analytic approaches are concerned, this update poses the greatest difficulties, 
since it produces the strongest correlations. In the case of the ASEP, this becomes obvious 
when looking at the model put on ring: All updates, except the parallel ]36|] an d particle- 



sequential JJTJ one, lead to a trivial state where correlations are absent. In the case of the 
parallel update, it is known |36|] that a particle-hole attraction appears: the probability 
to find an empty site in front of an occupied site is enhanced compared to the mean-field 
result. 

The model with parallel update and periodic boundary conditions has first been solved ex- 
actly in |8|,|36| using a cluster approximation (see below). Here the weights are decomposed 
into products of pairs of sites overlapping just one site. It turned out that this 2-cluster 
approximation becomes exact for periodic boundary conditions. 

We now make an analogous calculation for the model with open boundary conditions. The 
goal is to find the parameter set for which the density profile is flat. Note that a flat 
density profile does not necessarily mean that there are no correlations between the sites. 
We denote the probability for the pair configuration r i Tj +1 at site i, i + 1 by P Tin+1 (Tj = 1,0 
and i = 1, L). Let us assume that such a probability for a certain pair configuration is 
independent of the position of the pair. The condition P 01 = P 10 leads to a flat density 
profile ((ti) = (r 2 ) = ... = (tl}), but P TiTi+1 is not necessarily equal to P Ti ■ P n+1 - 
The 2-cluster approximation corresponds to a factorization of the weight P(ti, r 2 , tl): 

P(n, r 2 , t l ) = R T1 ■ P T1T2 ■ P T2TS • ... ■ P TL _ lTL ■ R TL , (51) 

where the Ks reflect the influence of the boundaries; they can be set to Rq = 1 and R± = r. 
It is sufficient to study a system of 3 sites. The generalization to larger systems is straight- 
forward |36| . 

The (exact) master equation for the stationary state reads x = Tx where x is a vector 
containing the 8 possible configurations, and T is the 8 x 8-dimensional transfer matrix for 
L = 3. By inserting Ansatz ([51]) into this master equation it is straightforward to show 
that (|5Tp is exact if and only if 

(l-a)(l-P) = l-p. (52) 

This means that the condition we have found for a constant density profile is exactly the 
same as for the other discrete-time updates, see (|5U[). The reason for this coincidence is 
still unkown. Note that (|52"D is not a simple mean-field line. 

In the following we will propose a MPA for the parallel update. ([51]) suggests that matrices 
denoting pairs of particles and/or holes should be used. The main problem is that the 
resulting algebra will be very complex, since it will contain many variables (matrices). On 
the other hand, if we use a MPA of the form fl2"7|) , a simple mean- field solution will not be 
found, e.g. there will be no scalar solution for the algebra. However, such an Ansatz could 
produce (|5T]) under condition (|52| ) in the form of a low- dimensional representation. 



Furthermore, even if a choice for the MPG is made, we still have to find a mechanism 
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which ensures stationarity. The transfer matrix, however, cannot be decomposed easily in 
products or sums of local terms. 

The main difference between the parallel update and the ordered- sequential update from 
the left to the right is that in the parallel update, a particle can move only one site to 
the right (per update of the chain). This enables us to use the ordered-sequential update, 
if we introduce a third state for particles that have been moved. The local (sequential) 
update operator then has to transform a third-state particle into a "normal" particle in 
the following update step. 
We write down the Ansatz 



\Po) 



(W\ 



5) 



which has to satisfy 



T\\ \P ) = \P ) 



(53) 



(54) 



The update-rule T for any pair of sites (i, i + 1) is now nine-dimensional. However, since 
the third state must not appear after the update, the last four rows and every third column 
are irrelevant (here set to zero). The explicit expression for T can be found in Appendix 
U We have 



R ■ Tl-i • ... • Ti • L 



(55) 



with 



L = £(g) 1(8) ...(g) 1, 
R = 1 <g) ... <g) 1 ®TZ, 
P = 1 ® 1... ® % <g> 1. 



(56) 



1. 



The mechanism for stationarity reads now 



T 




(w\c(d 





(57) 



*(f)|iO=(!)|v> 



with the new third-state matrix F. This leads to an algebra, which can be found in Ap- 
pendix 0. Note that the last bulk equation = DF excludes a scalar solution for the 
algebra. This is consistent with our earlier observation that there is no simple mean-field 
solution of the model. 

First, we can check the relations which connect the densities at the ends of the chain with 
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the current J. This calculation can be done for arbitrary system sizes and is presented in 
Appendix ^[ 

Second, it is possible to show that the algebra correctly describes a system of three sites. 
This can be done by taking the expressions for the weights of the eight possible configu- 
rations of the stationary state given by (Q) and by applying several times the algebraic 
rules. Thereby each weight can be expressed as a linear combination of the other weights^. 
The resulting system of linear equations turns out to be identical to those obtained from 
the transfer matrix T. 

Thirdly, it can be checked whether the algebra can be reduced to a generalized DEHP- 
algebra of the type 

pDE = aE + bD (58) 

with some numbers a, b. This equation induces certain relations between weights for a 
system of size L and L — 1, which can be checked using the exact solutions. It turned out 
that (BSD cannot be valid. 



We found a two-dimensional representation for the bulk algebra |38] , but despite intensive 
effort, we could not solve the complete algebra^. Therefore, Monte Carlo simulations have 
been performed p9|. The results will be presented in Section [12. 



4 Comparative study of physical quantities 

In the following we will investigate the consequences of the mathematical description de- 
veloped in the previous section. First we investigate the ASEP with periodic boundary 
conditions for the different updates. This will allow us later to distinguish between "pure" 
bulk effects and "boundary-induced" bulk effects. After that we will derive and compare 
the phase diagrams, density profiles and other physical quantities for the various updates. 



4.1 Periodic boundary conditions 

For random-sequential dynamics the stationary state of the ASEP with periodic boundary 
conditions is given by \Pq) = -^-Tr(A ® • • • ® A) where the elements E and D of the 
vector A satisfy the algebra (|30|) and the normalization is given by Zl = TrC L with 



C = E + D. Since the boundary equations (31) and fl32|) do not have to be considered, it 



is possible to find a one-dimensional representation of the matrices E and D which then 
become real numbers e, d. The current is calculated from = -^■ r Ti{C^~ l DEC L ^^ 1 ) = -, 
where we have used fl3~0|) and Zl = c L with c = e + d. In contrast to the case of open 



3 Note that formulas for macroscopic variables like the current are quite complicated for this small 
system; the current, for example, is a ratio of polynomials in a, and p containing 27 additive terms. 
4 This situation is very similar to |19[ where a MPA for the ASEP with a defect was proposed. 
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boundary conditions, the density p is now fixed and the density profile is constant, p = (tj), 
independent of j. With 



( Tj ) = ^Tr(C^DC L ^) (59) 

one finds p = d/c. Therefore the current is given by 

J(j>,p)=pp(l-p). (60) 

This is exactly the mean-field result, because in mean-field approximation a site is occupied 
with probability p and its right neighbour is empty with probability 1 — p. Since hopping 
then occurs with probability p, one obtains (|60|). 

This result is not surprising as the existence of a one-dimensional representation implies 
the absence of correlations between neighbouring sites, i.e. the MPA reduces to mean-field 
theory. 

We now turn to the ASEP with backward-sequential update on a ring. The state 

|P ) = -^TrlA®^ ® A] (61) 

is obviously translation- invariant and, because of (|42"D, stationary (the argument of the 
trace is a vector with a product of matrices in each component; the trace has to be applied 
to each component). The algebra can be directly solved with one- dimensional matrices 
(see Appendix [D]). The current J^(p,p) for the ordered-sequential update is not given by 



), but by| 



J_(j) = I-(W\C j - 1 DEC L - i - 1 \V) . (62) 
Zl 



For the case of periodic boundary conditions one findsf] 

J^(P,P)=PPT — — ■ (63) 
l-pp 

Fig. |3] illustrates this result. Again this result can be obtained directly using a mean- field 
argument. The site to the right of an occupied site is empty with probability 1 — (p— J<_). 
Here one has to take into account that the density is reduced by J<_ after the update 
of that site. Therefore the current satisfies J*_ = pp(l — (p — </<-)) which leads to (|63|). 
The maximal flow (p fixed) is reached for a density p™ ax (p) — -(1 — y/1 — p) > 1/2. The 
sequential update T<_ "likes" high p and high densities. The particle-hole symmetry can be 
used to determine these quantities for T_ simply by replacing p with 1 — p: "prefers" 

5 Note that eq. j6^ ) also applies for the sublattice-parallel update. 

6 The corresponding formula for the ASEP with hopping in both directions can be found in App. 



18 



high p and low densities. 
It is interesting to compare 



j) with the well-known result [36] for parallel update, 
J\\{p,p)=~-\ I 



l-App(l-p) . 



(64) 



For the parallel update a mean-field theory for the distances between consecutive particles 

becomes exact ||40|| . Similar results have been obtained recently by Evans ||37|| , who solved 

the ASEP with periodic boundary conditions in the presence of disorder^] with parallel 

update and a particle-sequential update by generalizing the approach of |TT]. 

J\\(p,p) is obviously symmetrical with respect to p = 1/2. It is maximal at pjj 1 ^ = 1/2 for 

all values of p, while p™ ax is always higher than 1/2 (except for p = 0). 

Hence, the maximal currents (for a given p) are J^(p max ,p) = ~(1 — ^1 — p) — 1 and 

J\\(p max ,p) = \{l - y/1 -p). It is intuitively clear that J^(p max ,p) > J||(p max ,p) holds, 
and it can be verified easily. 

Let us now return to the model with open boundary conditions. 



4.2 Phase diagram 

The phase diagram for random- sequential dynamics has been determined for p = 1 in fl3| ^ 
|T5| ]. This is no restriction since, as mentioned before, p only rescales time (and a — > a/p, 
(3 — > (3/p). Since it will turn out that the phase diagrams for the different updates are 
rather similar we will not repeat the results for the random-sequential update here. Instead 
we will first determine the phase diagram for the ASEP with open boundary conditions 



and discrete time by using the results of Section [D] and discuss the differences to the 
random- sequential case later. In Appendix |A| it is shown how the algebra ([P|) , fl4"5|) can be 
projected onto the DEHP-algebra. The representations of this algebra (and the resulting 
phase diagram) are known for all parameter values a, j3,p,q [^. . Therefore, we have 

obtained explicit expressions for the matrices and vectors, and have thus constructed the 
stationary state of our model; we do not have to be concerned about representations any 
more that might not satisfy ( |A.3| ) since the projection onto the DEHP-algebra already 
covers the whole parameter space. 

Furthermore, it is straightforward to calculate observables such as the current and the 
density, at least in principle. 

The mapping of App. [A| strongly suggests that the well-known phase diagram of the ASEP 
with random-sequential update and stochastic hopping in both directions |33|,|34|,|43| will 



also be valid for the ASEP with ordered- sequential update. This is indeed correct and has 



been proven directly in p9| . The phase diagram is shown in Fig. |j. The mean-field line 
is the curved dashed line. Also shown are density profiles calculated from Monte Carlo 
simulations. All well-known features of the phase diagram (high-density phase, low-density 



7 Each particle carries its own hopping probability pj, see also |4l|,p2| 



19 



phase, maximum current phase, coexistence line with linear density profile) are recovered. 
The intersection of the mean-field line (|50"D with the line a = (3 defines the endpoint a c 
(=Pc) of the coexistence line. This yields 




a c = 1 - VI - P- (65) 

In the case of deterministic hopping (p = 1), the maximum current phase vanishes, and we 
recover the result of Hinrichsen |26| . 

It is known that for the DEHP-algebra, the mean-field expressions for the current are 
exact. Since the mean-field line touches all three phases, we can calculate the corresponding 
currents and bulk densities for our model. Our results are 

in the low-density phase, 
J{a,f3,p) = ^ ffzf m the high-density phase, (66) 

in the maximum current phase, 

which is in excellent agreement with our numerical data. 

Since the relation J(a,(3,p) = flp(a, 0,p,x = L) is exact, we immediately get the bulk 
density in the high-density phase: ptmik = pfzf- The bulk density in the low-density phase 
will be obtained below. Fig. ^| shows a space-time diagram for a point on the coexistence 
line (p = 0.75, a = (3 = 0.3) produced by the Monte Carlo simulation. The well-known 
fluctuating shock can nicely be observed. 

Typical space-time diagrams for the different phases can be found in Fig. ^. The "jams" 
in the high-density phase move from the right to the left. In the low-density phase, groups 
of particles (small jams) move from the left to the right. 

For the case of parallel dynamics extensive Monte Carlo simulations have been performed 
3~9fl which revealed that the phase diagram looks essentially the same as before (Fig. |j). 



Again, we find a high- (low-) density and a maximum current phase and a linear density 
profile along the line a = (3, until the 2-cluster line fl5"2"D is touched. 

Since we have a particle-hole symmetry in the model, the density for a = (3 and odd lattice 
sizes has to be 1/2 for site . Therefore, the 'bulk density' p(a,0,p,x = ) in the 
maximum current phase must be 1/2. 

Equations fl5"TD and fl52"D can be used to calculate the current^ 

a ^E^2 m the low-density phase, 
J(a,j3,p) = < f3^z§2 in the high-density phase, (67) 



v-P 2 



in the maximum current phase, 



and the remaining bulk densities ptmik = ^z^r- (low-density) and p bu ik = j^E§2 (high- 
density phase). These results are exact on the 2-cluster line fl52"D. For general values of 
a and (3 we find an excellent agreement with the numerical data. Fig. ^| presents typical 
space-time diagrams in the various phases. The particle-hole attraction is apparent. 

8 For the parallel update the current is given by (|3 
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Note that for every update the maximal current J(p max ,p) on the ring (as calculated in 
Section |4.1|) is equal to the value of the current in the maximum current phase for that 
update. The value of p max gives the corresponding bulk density. This shows that the 
current in the maximum current phase is determined by the 'capacity' of the ring. 



4.3 Density profiles and phase transitions 

For the parallel update we performed extensive Monte Carlo simulations in order to de- 
termine quantities like density profiles and correlation functions. In this section we will 
compare these results with those obtained analytically and numerically for the other update 
types. 



The results presented in Section [4.2| show that one can distinguish at least three different 
phases with respect to the flow. In the low- (high-) density phase the current depends 
for a given value of p only on a (/3) and in the maximum current phase the current is 
independent of the chosen in- and output rates. This basic structure of the phase diagram 
is the same for all types of update although the bulk properties change drastically. 



A more detailed analysis of ASEP with random-sequential update has shown |TJ| that 



the system is governed by two independent length scales £ Q and which represent the 
influence of the boundaries. At the critical value of a and j3 these length scales diverge 
and a phase transition occurs. This divergence produces two additional phases compared 
to the mean-field results. Thus the parameter space is divided into five different phases. 
The low-density phases AI and All, where the bulk properties are determined by the value 
of a, the high-density phases BI and BII, where the output rate f3 determines the flow, 
and the maximum current phase where the flow is independent of a and (3. The phases 
AI and All (BI and BII) are distinguished by the behaviour of the density profile near the 
boundaries (see below). 

We checked this scenario for the parallel and ordered-sequential update. In order to avoid 
difficulties due to the p-dependent scale factor of the phase diagram we analyzed the density 
profiles for p = 0.75. Our simulation results are obtained for systems with 320 sites. Since 
finite-size corrections are rather small, this size is already sufficient to obtain the behaviour 
in the limit of large L. 



We first measured the density profile in the low-density phases AI and AIL In |14| , [L5[ it has 
been shown that for large L the asymptotic behaviour near the boundaries changes from 
a pure exponential decay in the phase AI to the enhanced exponential decay according 
to exp(— (L — x)/^a)/\/L — x in the phase All. This change of the behaviour near the 
boundaries is due to the divergence of the length scale £^ at the transition line. £p remains 
infinite throughout the whole phase AIL Unfortunately, the divergence of cannot be 
calculated directly because the relevant length scale £, which determines the exponential 
decay in AI, is given by 

r'Hc'-G 1 !- ( 68 ) 
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Thus the correlation length £ remains finite at the transition line from AI to AIL In 
contrast to the divergence of £p, the asymptotics of the density profile can be checked 
against the numerical results. Fig. ^ shows that the density profile for a = 0.40, (3 = 
0.42 and parallel update can be nicely fitted by a pure exponential decay, but within the 
phase All the enhanced exponential form has to be used (see also the insert in Fig. 
Another characteristic line crossing the low-density phase All is the mean-field line, where 
the density profile is completely flat. This line separates the monotonously decreasing 
density profiles from monotonously increasing profiles, but the asymptotic behaviour is 
left unchanged. The behaviour of the density profiles in the high-density phases BI and 
BII can be obtained using the particle symmetry of the model. 

The transition between the high- and low-density phases is driven by the diverging length 
scale £. This divergence occurs although both length scales are finite, because at the 
transition line £ a and £p coincide. The qualitative agreement of the density profile strongly 
suggest that the relevant length scale in the low-density phase is given by fl68|) also for 
discrete-time updates. Moreover, the comparison of the numerically estimated correlation 
length shows that the correlation lengths of the parallel and ordered-sequential update are 
identical and differ from those obtained for the random-sequential only by a constant factor 
(Fig. |8|). Exactly at the transition line a = (3, one finds a linear density profile as shown 
in Fig. |]. This linear profile is result of a fluctuating shock front which separates for a 
given time a high-density region from a low-density region. The position of the shock front 
fluctuates through the whole lattice such that a time average over all single time profiles 
gives a linear profile. 

The transition from the low- (high-) density phase All (BII) to the maximum current phase 
is characterized by a change of the asymptotic behaviour of the density profile from the 
enhanced exponential to a pure algebraic decay. The transition is driven by the divergence 
of the length scale £ a (£g) for a — > a c ((3 — > /3 C ). Both length scales are infinite not only 
at the transition line but for all values of a and /3 larger than the critical value. Therefore 
an algebraic decay of the density profile can be observed in the whole maximum current 
phase. Fig. [TI] shows the diverging correlation length for different update types. Again the 
length scales of the updates with discrete time agree and the length scales produced by 
the random- sequential update are larger if one considers the same in- and output rates. 

The flow is generated by the (lO)-clusters (the mobile pairs), while the other 2-cluster 
configurations exclude hopping of particles. Therefore, we measured density profiles of the 
probabilities P Tin+1 of 2-cluster configurations 7V7i+i. One gets a flat profile of the mobile 
pairs PiiQ i+1 (see Fig. [□]) which is a consequence of eq. Q5BD for the current in the case of 
parallel dynamics^. The nontrivial part of the density profile is produced by the immobile 
pairs. The identity Pi l o i+1 = PoAi+D which is already known from the periodic system ||, 
is only true for flat density profiles. One observes qualitatively the same behaviour for the 
random-sequential, but not for the ordered-sequential update. 

9 This is not true for the ordered-sequential update, since in this case the current depends on the local 
defect generated by the update. Hence, the current depends on DE (see (|62|)) instead of DE (see (]33|)), 
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For the ordered- sequential update the transport of a local defect changes the behaviour of 
the P TiTi+1 profiles: The density of particles in front of hole at the defect site (in the case 
of backward-sequential update) determines the flow and therefore none of the 2-cluster 
configurations is in general translation-invariant after a sweep through the lattice. 



5 Discussion 

We have presented an extensive comparative study of the ASEP with different updates. 
The purpose of this investigation is threefold: First of all, despite the importance of the 
ASEP, not much has been known about its properties for discrete-time updates. Second, 
we tried to obtain a better understanding of the similarities and differences of the different 
updates. Finally, there is also a practical aspect. In |?5],f0]] it has been shown that 



different dynamics perform quite differently in Monte Carlo simulations. In order to save 
computational time it might therefore be useful to choose a certain update. Then it is also 
necessary to know how to translate the results into those for other updating schemes. 

The main tool that we used in our investigations was the MPA which allows for an analyt- 
ical solution for the cases of random-sequential, ordered-sequential and sublattice-parallel 
updates. We also proposed a MPA for the important case of parallel dynamics, but unfor- 
tunately we were not able to find a general representation of the resulting matrix algebra. 
Therefore extensive Monte Carlo simulations have been performed in order to determine 
the phase diagram. These numerical results, together with an analytical solution for a spe- 
cial line in parameter space, allowed us to conjecture analytical expressions for the phase 
diagram. 

Our results show that the phase diagram has the same basic structure for all the updates 
investigated here. One finds three different phases characterized by the value of the current. 
For a > (3 and (3 < /3 c (p) the system is in the so-called high-density phase. Here the current 
depends only on the removal rate j3 since particles are inserted much more efficiently than 
they are removed. Just the opposite situation is found in the low-density phase, a < (3 
and a < a c (p). Here the removal is much more effective than the insertion and the current 
depends only on a. Note that for each update a c (p) = f3 c {p) and the functional forms of 
the currents in the high- and low-density phases are identical (see Table 1). 

Finally, for a > a c (p) and (3 > (3 c (p) one finds the maximum current phase. Here the 
current is independent of a and (3. Both insertion and removal are so effective that the 
current is only limited by the "bulk capacity" . Indeed, the current in this phase is identical 
to the maximal current of the corresponding system with periodic boundary conditions. 

Phase transitions between these phases are driven by diverging length scales or £p 
within the high- and low-density phases. These length scales depend on the rates a and 
(3, respectively. In contrast to that, the periodic systems exhibit only extremely short- 
ranged correlations. The strongest correlations are found for the parallel update, but even 
here already the 2-cluster approximation is exact. This leads to exponentially decaying 
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correlation functions with a rather short correlation length (except for p « 1/2 and p ~ 1 ) 



45| . So the long-ranged correlations found for the open system are due to the boundary 
conditions and the transitions are genuine "boundary-induced phase transitions". This 
makes it also understandable why the phase diagrams for the different updates look so 
similar, the only difference being the location of the transitions and the functional form 
of the observables. Furthermore the numerical analysis shows that the update type does 
not change the "universal" properties of the model: We observe the same asymptotics of 
the density profiles in the different phases and also the qualitative behaviour of the length 
scales near the phase boundaries does not depend on the update. For the discrete updates 
the length scales agree even quantitatively. 

Recently, for the case of periodic boundary conditions there has been some progress in 
the understanding of the differences between random-sequential and parallel update p6 



For the latter so-called "Garden of Eden" (GoE) states^ exist. These states can not be 
reached dynamically, they do not have a predecessor. By eliminating these states one finds 
in the reduced configuration space that mean- field theory becomes exact, as it is for the 
case of random-sequential update (but here in the full configuration space since Garden 
of Eden states do not exist). Therefore the existence of GoE states is the reason for the 
correlations in the ASEP with parallel dynamics and periodic boundary conditions. Since 
the bulk dynamics for the ASEP with open boundaries is the same as for the periodic case, 
we expect the GoE states to play an important role also in that case. 

In this paper, valuable information has been gained by means of low- dimensional repre- 
sentations of the matrix algebras. A one-dimensional representation clearly corresponds to 
a simple mean-field approximation. However, in the case of two (or higher) dimensional 
representations, nothing is known about the physical interpretation of these solutions. For 
example, it could be possible that there is a close connection between cluster approxi- 
mations and these representations; we remark that it is possible to write the expectation 
values of densities and correlations of any exact (stationary) solution of a 2-cluster approx- 
imation as a product of two-dimensional matrices precisely in the form resulting from a 
MPA. 

Recent investigations have shown that the ASEP is capable of reproducing the essential 
features of traffic in a city such as Geneva [[|7]]. The authors studied a simple extension of 
the one-dimensional ASEP to two dimensions^} Therefore, it would be most interesting to 
generalize the MPA to higher dimensions^. A first step would be to study the ASEP on 
a ladder. The analytical method used in this paper is directly applicable to "stochastic" 
ladders j52J . Since the groundstates of certain quantum spin ladders have been constructed 



recently J53],^4] using optimum groundstates (see Section |2.3| ), there seems to be a good 
chance to find low-dimensional representations of the corresponding algebras. This would 



10 Here on has to distinguish between particles which moved in the previous timestep (velocity 1) and 
particle that did not move (velocity 0). 

11 In (4^ |5^] other generalizations of the ASEP have been used for simulations of urban traffic in Duisburg 
and Dallas. 

12 See |5ll] for a numerical investigation of a two-dimensional ASEP. 
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lead to a better analytical understanding of the fascinating phenomena that occur in higher- 
dimensional nonequilibrium systems |2j,]55|. 
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A Mapping of the ordered-sequential algebra onto the 
DEHP-algebra 

We here treat the general case where particles are also allowed to hop to empty sites on 
the left with probability q. In this case the bulk algebra (|44] ) is generalized to 

[E,E) = [D,D] = 0, 

(1 — q)ED + pDE = ED, (A.l) 

qED + {l-p)DE = DE, 

with the boundary conditions (flop . We first note that 

[E + D, E + D] = (A.2) 

holds for all values of p, q. 

The crucial step is now to demand 

E = E + Xl, 

D = D-Xl (A.3) 

with some (real) number A. Note that this is the simplest way to satisfy ( |A.2| ). 

Now, we must show that a solution (representation) of the algebra (^4|), fl45[) plus equations 

( |A.3j ) exists for all possible values of the parameters a,j3,p,q. 

( |A.3j ) reduces the algebra from eight equations to three equations: 

pDE-qED = X(l-q)E + X(l-p)D, 

a(W\E = \(l-a){W\, (A.4) 
0D\V) = X\V). 



We define 



D := X(l-p)D, 

E := X(l-q)E, (A.5) 



A 2 



(l-g)(l-p) 
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and rewrite (|A.4 ) as 



pDE - qED 
a(l-p)(W\E 
(3(l-q)D\V) 



E + D, 

(l-a){W\ 

\V). 



(A.6) 



This is the algebra for the ASEP with random- sequential update and hopping in both 
directions (with probability p and q, respectively), but with the same local transfer matrix 
and the same boundary conditions as in our model. The algebra was solved for p = 1 and 
q = by Derrida et al. (T^j with infinite-dimensional matrices. Note that the vectors {W\ 
and \V) of their solution have to be rescaled with and respectively. A thorough 



discussion of the algebra ( |A.6|) can be found in |33|, |34j,|j|,|56|. 

We write down an explicit representation of the algebra (|A.4j ) with A 



D 



1 

P 



( f a x 

110 

11 

1 



\ 



1 in the case q = 0: 



(A.7) 



E 



1 

P 



/ p(l-a) 
' a 

a 2 







1 — p 
1 — p 







1 — p 
1 — p 







1 — p 



(A. 



(W\ = (1,0,0,0,...), \V) 



(l\ 





VJ 



(A.9) 



o 1 a 2 = ^-[(l-p)-(l-a)(l- / 9)]. 



(A.10) 



As expected, constraint ( p0| ) leads to an effectively one-dimensional representation. 
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B Symmetries of the density profiles for the ordered- 
sequential updates 



As shown in Section y the stationary state for T_> is simply given by (|4q) . When inserting 
( |A.5| ) into (^6|), we get a connection between the density profiles^ of the stationary states 
produced by and T_>: 

p^(a,f3,p,x) = p^(a,P,p,x) - A— — . (B.l) 

One also has J(a,j3,p) = j3p(a,j3,x = L) which leads to 

p-+(a,p,p,x) = p^(a,P,p,x) -J(a,P,p) . (B.2) 

(a, /3, p, x) is therefore always lower than p^_ (a, (3, p, x). The current J is not x-dependent. 
This means that the density profile of the ordered-sequential model for a given set of pa- 
rameters a, /3,p is, up to a constant (the current), the same for both directions of the order 
of update. The stationary states produced by T_» and T<_ for a given set of parameters are 
always in the same phase. It is obvious that the "crucial step" ( |A.3|) is the simplest way 
to obtain such a behaviour. 
When using the particle-hole symmetry 

P-*{(x,f3,p,x) = 1 - p^.{(3,a,p,L - x + 1) , (B.3) 

we immediately see that the density profile on the line a = f3 has to be symmetric with 
respect to point L/2. For this case, we further obtain p^_(a,p, = 1+ J ( a ' p -> . Finally, by 
using ( |B.2j ) and ( |B.3j ), and the results for the current and ptmik m the high-density region, 
the bulk density in the low-density region ptmik = a/p is obtained. 



C Symmetric diffusion 

We briefly discuss the case of symmetric diffusion [p = q) for the ordered-sequential update. 
The density p(a, j3,p,x) at site x (1 < x < L) is given by 

p(a,P,p,x) = -^(WIC^DC^IV). (C.l) 



flOl) and flOp yield 



[D,C\ = [D,E + D] = -(E + D) = ^C. (C.2) 



13 In the following we will denote the local density (t x ) at site x by p(a,(3,p,x) in order to stress the 
dependence on the other parameters. 
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This means that the density can be calculated immediately by commuting D with the 
help of (|C.2|) through the chain in (|C.1|) . By using the boundary conditions, the density 



can easily be calculated, which in turn makes it possible to estimate the current [29]. It 
is intuitively clear that the current vanishes in the thermodynamical limit L — > oo for 
arbitrary values of a,f3,p. 



D ASEP with ordered-sequential update on a ring 

Since we expect that the periodic system can be described by an one-dimensional repre- 
sentation of the algebra ( |A.1[ ), we are looking for solutions of ( |A.1| ) with real numbers 
6 . d ^ 6j dt 



(1 — q)ed + pde = ed , (D.l) 
qed+(l-p)de = de . (D.2) 



The normalized density p is given by 

d{e + d) L - 2 {e + d) 1 



{e + d) L - l {e + d) (e + d)' 



d- -. (D.3) 



There is some freedom how to choose e, d, e, d. We set e + d — e + d — 1. Therefore, ( D.l ) 
yields 

, = P fa-i) + i- 9 

p{q-p) + l-q 
The current J is given by J = pde — qed and one obtains 

p(q-p) + l-q 

Obviously, J(p, q) = for the case of symmetric diffusion, p = q. 

Finally, we like to point out that the fact that the periodic system is described by a 
one-dimensional representation implies that in this case mean-field theory is exact. 



E T]i in a three-state notation 

The operators 7Z,£ can be written as: 

/ 1 (3 

ft = 1-/3 1 I , £=| 1 I . (E.l) 
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The basis chosen for 1Z and £ is (0, 1, —1). 

The local transfer-matrix % for any pair of sites + 1) is nine-dimensional. Since the 
third state (—1) may not appear after the update of the whole chain, the last four rows 
and every third column are irrelevant and here set to zero: 



/ 1 























\ 





1 
































V 


























1 — p 








1 




















1 








1 













































































































0/ 



(E.2) 



The basis is (00, 01, - 1, 10, 11, 1 - 1, -10, -11,-1-1). 



F Matrix algebra for parallel dynamics 



The MPA for parallel update proposed in Section 3.2 leads to the following algebra 



and the boundary conditions 



EE 
ED 
pt)E 

[1 — p)DE + EE 
DD + FD 




EE, 
ED, 
EE, 
DE, 
DD, 
DF, 



(W\E(l-a) = (W\E 
(W\D = (W\D 
(W\aE = (W\F 
E + (3D\V) = E\V) 
[l-P)b + F\V) = D\V) 



[FA) 

;f.2) 
;f.s) 
;f.4) 
;f.s) 
;f.6) 



(F.7) 
(F.8) 
(F.9) 
(F.10) 
(F.ll) 



Note that the last bulk equation = DF excludes a scalar solution for the algebra. This 
is consistent with our earlier observation that there is no simple mean-field solution of the 
model. 



29 



G Check of current-density relations 

Defining C = E + D and C = E + D + E it is easy to see that 

[C,C] = (G.l) 

holds. This is implied by probability conservation (the columns of % add up to one) and 
the exchange mechanism fl57|). The boundary equations lead to 



(W\C = (W\C (G.2) 

and 

C\V) = C\V) . (G.3) 

The bulk equations flFJ), flFJ) give 

CD = ED + DD = DD + ED + ED = CD . (G.4) 

We can now check the relations which connect the densities at the ends of the chain and 
the current J which has to be constant throughout the chain: 

J(a,j3,p) = j3p(a, j3,p,x = L) (G.5) 
J(a,/3,p) = a[l — p(a, /3,p, x — 1)] . (G.6) 

Note that the first (second) equation would not hold for the sequential update from the 
left (right) to the right (left) and that in fact we have only to prove one of these equations 
because we can make use of the particle-hole symmetry of the model. We therefore write, 
using the algebra (equations (|F.11|) , ([F . 3|) ) : 

p(a,t3,p,z = L) = ^-(W\C L - l D\V) (G.7) 
= f3)(W\C L - l D\V) +^ r p(W\C L - 2 DE\V). 



Making use of ( |G.4|) , commuting C to the left end of the chain (|G . 1| ) , transforming it to 
C there Q, and using J(a,/3,p) = ^p(W\C L - 2 DE\V) we get 

p(a,(3,p,x = L) = -L(1-/?)<W|C^D|V0 +.L P (W\C L - 2 DE\V) 

= (l-(3)p(a,P,p,x = L) + J(a,f3,p), (G.8) 

which is the desired result. 
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Table 



rand. -sequential ordered-seq. (T<_) parallel 



low-density 
pnase 


J 


= pa{l — a) 

p = Oi. 


J = 
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a p—a 
p 1— a 
- ^ 


J = a p - a 2 

p—a z 
_ a(l-a) 
• p—a 2 


high- density 
phase 
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= P p(i-P) 
p=i-P 


J = 
P = 


. /3P-/3 
■ p 1-/3 
. 1P-/3 
pl-/3 


J = P V ~L 
o 

P ~ p-p 


max. current 
phase 




J = l 
p = \ 


J = 

PL/2 = 


i- v / r r p 
i 


T _ i--/i=p 

PL/2 = | 














critical rate 




UL C 2 


a c — 1 


- y/l-p 





Table 1: Comparison of currents and (bulk) densities in the three phases for the different 
updates. The bulk densities for T_> are given by p — J for T<_. The currents for and 
TU are identical. 
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Figure Captions 



Fig. 1 Definition of the ASEP. 

Fig. 2 (i) Ordered-sequential update from the right to the left, (ii) from the left to the 
right, (iii) and the sublattice-parallel update. 

Fig. 3 Fundamental diagram (flow vs. density) for the ASEP with ordered-sequential up- 
date and periodic boundary conditions. The vertical lines indicate the location of 

p™ ax (jo). 

Fig. 4 Phase diagram for the ASEP with ordered-sequential update TL for p = 0.5. 
The mean-field line (|50|) is curved dashed line. The straight dashed lines are the 
boundaries between the phase AI and All (BI and BII). The inserts show density 
profiles obtained from Monte Carlo simulations. 

Fig. 5 The space-time diagram for p = 0.75, a = (3 = 0.3. The microscopic shock moves 
freely on the chain. 

Fig. 6 Typical space-time diagrams for the various phases. 

Fig. 7 Density profiles for parallel update (p=0.75) in the low-density phases AI and AIL 
The insert compares the exact asymptotic form with a pure exponential decay. 

Fig. 8 Log-Log plot of the correlation length for different types of updates (OS=ordered- 
sequential, PARA=parallel, RS=random-sequential). The length scales are obtained 
from an exponential fit of the density profile, a and (3 are chosen such that a+/3 = 0.6 
holds (a c = 0.3). 

Fig. 9 Density profiles near the first order transition at a = (3 = 0.3 for the random 
sequential update using the exact results of Again a + (3 = 0.6 holds. 

Fig. 10 Divergence of the correlation length £ Q near the transition to the maximum current 
phase (J3 = 0.75). 

Fig. 11 Density profile of the pair probabilities P rm+1 for the parallel update in the low 
density phase All (a = 0.40,/? = 0.75, p = 0.75). 
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Figure 1: Definition of the ASEP. 
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(ii) 



(iii) 






V_J< V_J< V_J< V_J< 



Figure 2: (i) Ordered-sequential update from the right to the left, (ii) from the left to the 
right, (iii) and the sublattice-parallel update. 
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Figure 3: Fundamental diagram (flow vs. density) for the ASEP with ordered- sequential up- 
date and periodic boundary conditions. The vertical lines indicate the location of p™ ax (p). 
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Figure 4: Phase diagram for the ASEP with ordered-sequential update T<_ for p = 0.5. The 
mean-field line ( |50l) is curved dashed line. The straight dashed lines are the boundaries 
between the phase AI and All (BI and BII). The inserts show density profiles obtained 
from Monte Carlo simulations. 
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■ moving particles 
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Figure 5: The space-time diagram for p = 0.75, a — (3 — 0.3. The microscopic shock moves 
freely on the chain. 
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Figure 6: Typical space-time diagrams for the various phases. 




Figure 7: Density profiles for parallel update (p=0.75) in the low-density phases AI and 
AIL The insert compares the exact asymptotic form with a pure exponential decay. 
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Figure 8: Log-Log plot of the correlation length for different types of updates (OS=ordered- 
sequential, PARA=parallel, RS=random-sequential). The length scales are obtained from 
an exponential fit of the density profile, a and (5 are chosen such that a + [3 = 0.6 holds 
(a c = 0.3). 
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Figure 10: Divergence of the correlation length £ Q near the transition to the maximum 
current phase {(5 = 0.75). 
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Figure 11: Density profile of the pair probabilities P n r i+1 for the parallel update in the low 
density phase All (a = 0.40, (5 = 0.75, p = 0.75). 
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